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RESULTS  ON  THE  DETECTION  OF  SIGNALS  IN 
SPHERICALLY  INVARIANT  RANDOM  NOISE 


INTRODUCTION 

Vershik  [1]  introduced  the  class  of  spherically  invariant  random  processes  (SIRP)  as  part  of  a  study  of 
properties  exhibited  by  Gaussian  random  processes.  Since  that  introduction,  papers  exploring  various 
aspects  of  this  class  of  processes  have  appeared.  These  studies  have  examined  questions  regarding 
representations  of  such  processes  [2-5]  as  well  as  estimation  and  detection  in  such  processes  [6-8], 
Recently,  this  work  has  been  applied  to  the  problem  of  modeling  certain  types  radar  clutter  processes  [9, 10] . 

In  this  study,  we  reexamine  the  structure  of  the  optimal  processor,  i.e.,  the  likelihood  ratio,  for 
detecting  a  known  signal  in  additive  noise  modeled  as  a  SIRP.  The  structure  of  this  optimal  detector  is  that  of 
the  estimator-correlator,  a  result  that  follows  by  direct  application  of  a  result  attributable  to  Schwartz  [11], 
However,  for  the  sake  of  clarity,  a  simple  derivation  of  this  result  is  presented.  In  addition,  we  discuss  how 
the  estimator-correlator  structure  leads  to  suboptimal  structures  that  have  been  obtained  by  previous 
investigators.  We  also  briefly  point  out  that  the  locally  optimum  detector  for  this  problem  consists  of  the 
same  optimal  estimator  used  in  a  different  manner  than  in  the  optimal  detector.  Finally,  we  derive  an 
alternate  formulation  of  the  likelihood  ratio  and  the  optimal  estimator  that  appears  in  the 
estimator-correlator  structure.  This  alternate  formulation  is  important  because  it  allows  the  optimal  estimate 
to  be  computed  without  knowledge  of  a  prior  density  for  the  random  parameter.  Instead,  knowledge  of  a 
certain  marginal  probability  density  function  (pdf)  is  required.  In  practice,  this  pdf  is  more  likely  to  be 
available  to  the  investigator  than  is  the  prior  density. 

STRUCTURE  OF  THE  OPTIMAL  DETECTOR 

In  this  section,  the  optimal  detection  structure  for  detecting  an  additive  signal  in  a  spherically  invariant 
random  process  is  derived  and  shown  to  be  a  general  structure  applicable  to  the  whole  class  of  processes. 
The  resulting  structure  is  a  function  of  an  optimal  estimator  of  a  random  quantity.  This  structure  reveals  an 
intimate  relationship  between  optimal  detection  and  optimal  estimation  for  this  class  of  processes. 

Consider  the  following  hypothesis  test: 

H0:  x  =  y 


Hi',  x  =  y  +  s 


where 

y  is  an  m-dimensional  complex  noise  vector 
s  is  an  m-dimensional  complex  noise  complex  vector 

=  (is,  (3  =  yeJ(t>  =  complex  number. 
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From  this  point  only  complex  processes  are  examined.  However  the  corresponding  results  for  real 
processes  are  easily  obtainable  from  our  approach. 

If  a:  is  a  sample  from  a  SIRP,  then  x  has  the  following  pdf  under  hypothesis  i  (i  =  0,1): 

fix)  =  C^j*|  Jo"  7T  exp  {’  S)  fMdr'  <» 

where 

qt  =  (. x-s y  (x-s),  i  =  1 

=  i  =  0 


$  is  the  normalized  correlation  matrix 
r  is  the  variance  of  underlying  Gaussian  process 
fff)  =  pdf  of  r 
1 1  is  the  matrix  determinant 
t  is  the  complex  conjugate  transpose. 


Clearly  Eq.  (1)  represents  the  pdf  of  a  conditional  Gaussian  process  with  random  variance  t,  averaged  over 
the  variation  in  t.  Close  examination  of  the  conditional  Gaussian  pdf  shows  that  it  is  a  multivariate  member 
of  the  exponential  family,  i.e.,  the  pdf  is  of  the  form 


where  in  our  case 


fi(x\r)  =  c(r)h(x)exp 


i 

-  ^2  V j(r)Tj(x)  >, 
j-  1 


(2) 


c(t)  -  - 

(27T)m  |^>|/" 

h(x)  =  1 
/  =  1 

Vi  (?)  =  — 

T 

TiCx)  is  a  sufficient  statistic  for 

Our  first  result  follows  by  direct  application  of  a  general  result  given  by  Schwartz  [11],  However,  we  now 
give  a  brief  derivation. 
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If  we  reparameterize  our  conditional  Gaussian  pdf  by  setting  a  =  1/r ,  we  obtain 


(3) 


where  we  have  substituted  Zi  =  qj2 .  If  we  now  differentiate  Eq.  (3)  by  z,-,  we  obtain 


(4) 


Since  the  minimum  mean-square  estimate  of  a  is  given  by  the  conditional  mean  estimate  (CME)  [12], 


(5) 


and  we  immediately  obtain 


(6a) 


dlnfi 


(6b) 


At  this  point  we  note  that  we  may  use  a  dummy  variable  z  in  place  of  z,  in  our  formulation,  i.e. ,  we  now  use 
Z  =  qi/2  instead  of  Z;  =  qt/2  .  From  Eq.  6(b)  we  then  obtain 


(7a) 


where 


(7a) 


Since  under  either  hypothesis  we  have^x)  =  g(qt),  i.e. ,  the  subscript  i  indicates  that  we  have  a  function  of 
qt ,  and  since  we  are  using  z  =  qt/2,  it  is  easy  to  show  that  the  form  of  £f[a|z]  is  the  same  under  both 
hypotheses.  We  denote  this  form  E[a\z]  and  differentiate  between  the  two  hypotheses  by  evaluating  this 
form  at  the  different  points,  z  =  qi/2 ,i  =  0, 1.  Therefore,  since  -ln^]z.=0  is  the  same  under  both 
hypotheses,  we  obtain  finally 


(8) 


The  structure  of  the  likelihood  ratio  is  therefore  given  by 


(9) 


Equation  (9)  shows  explicitly  the  relationship  between  optimal  estimation  and  optimal  detection  for  this  class 


of  processes. 
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A  GEOMETRIC/FUNCTIONAL  INTERPRETATION  OF  THE 
OPTIMAL  DETECTION  STRUCTURE 

In  this  section,  we  interpret  the  optimal  detection  structure  as  the  area  under  a  curve  in  a  certain 
coordinate  system.  This  interpretation  then  leads  to  the  conclusion  that  the  optimal  detection  structure  is 
equivalent  to  a  comparison  of  the  matched  filter  with  a  variable  threshold,  which  is  a  function  of  qo ,  the 
quadratic  form  in  x,  and  the  detection  threshold  T that  yields  the  desired  probability  of  false  alarm. 

From  Eq.  (8),  we  may  write  the  optimal  detection  structure  as 

Mon 

U/2E[a^z]dz  <  T •  (10) 

Ho 

Note  that  since  is  in  general  a  positive  definite  matrix,  the  limits  of  integration  are  always  positive. 
Let  us  assume  temporarily  that  q q  >  q\  and  T  >  0.  Also,  since  a  is  a  nonnegative  random  vari¬ 
able,  the  integrand  E[a\z],  which  is  a  function  of  z,  is  also  nonnegative  for  positive  z.  Thus,  the 
structure  in  Eq.  (10)  compares  the  area  under  a  curve  in  the  E[a  \  z]  —  z  coordinate  system  to  a 
threshold  T.  Figure  1  shows  this  comparison. 


N 

X 

LU 


q  /  2  q  /  2 

1  0 

Fig.  1  —  Geometric  interpretation  of  optimal  detector 

Now,  since  £[g!|z]  is  nonnegative  for  positive  z  and  the  limits  q\/2  and  qo/2  are  always  positive,  this 
area  is  <  T  if  and  only  if  q\/2  >  qo/2  -  K,  i.e.,  q\/2  is  closer  than  K  units  to  qo/2.  Note,  however,  that 
unless  E[o!  |e]  is  a  constant  for  positive  z  the  number  of  units  K  is  actually  a  function  of  the  location  of  qo/2  on 
the  z  axis.  Thus,  the  comparison  indicated  in  Eq.  (10)  is  equivalent  to  the  comparison 

Hx 

Go  q\  >  r,  rj^ 

2  2  <  f ^ 0,T),  (11) 

Ho 

where  f(qo,  T)  is  a  variable  threshold  that  is  a  function  of  both  qQ  and  T.  After  substitution  for  the 
quadratic  forms,  Eq.  (11)  becomes 


AREA  COMPARED  TO  T 
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Re(s‘*-lx)  >  flq0,  7)  +  |  s'&'s,  (12) 

Ho 

which  is  a  matched  filter  compared  to  a  variable  threshold.  Theoretically  this  function  f{qo,  7)  may  be 
obtained  by  finding  the  function  f(qo)  that  satisfies 

f  <?i  /2 

l  E[a\z]dz=T.  nT. 

J  <?o/2-/('?o)  (13) 

At  this  point,  however,  we  abandon  the  geometric  interpretation  of  this  detection  structure  and  simply 
examine  the  structure  given  by  Eq.  (10)  from  a  functional  point  of  view.  Define  a  function  F(q)  as 

A 

F (<?)  =  r  E[a\z]dz.  (14) 

Since,  as  we  have  previously  observed,  E[o;|x]  is  nonnegative  for  positive  z  and  q  is  positive,  the  function 
F(q)  that  we  have  defined  is  a  monotonically  increasing  function  of  q.  In  this  new  notation,  we  rewrite  the 
log-likelihood  test  given  in  Eq.  (10)  as 


or  equivalently, 


Hi 

F(.qo)-F(<h)  <  T, 
Ho 


Hi 

F(q0)-T  >  F(q{). 
Ho 

Nov  ,  since  F(q)  is  monotonically  increasing,  we  obtain 


Finally,  since 


Hi 

r\F(q0)-T)  >  qx. 

Ho 

qx  =  x'^'x  -  2  Refs'^x)  + 


=  qo~  2  Re(^'<I>  *x)  + 

we  obtain  after  substitution  in  Eq.  (16) 


Hi 

Re(s'$_1x)  ^ 

Ho 


f-\F\F(q0)-T)+±s'*-'s. 


(15a) 


(15b) 


(16) 


(17) 


(18) 
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Comparison  of  this  result  with  the  result  in  Eq.  (12)  reveals  that 

Aqo,T)=^-±F-\F(q0)-D.  (19) 

This  equation  may  or  may  not  be  solvable  for  a  specific  problem.  However,  the  main  advantage  to  this 
interpretation  is  the  insight  it  gives  into  the  operation  of  the  optimal  detector.  Note  that  in  this  functional 
interpretation,  we  made  no  assumptions  about  either  the  ordering  of  q0  and  q\  or  the  sign  of  T. 

Observe  that  if,  in  Eq.  (10),  T  =  0,  which  corresponds  to  a  unity  threshold  for  the  likelihood  ratio,  then 
the  log-likelihood  ratio  reduces  to 


> 

<7o  <  qu  (20) 

H0 

which  in  turn  is  equivalent  to  the  matched  filter  and  is  independent  of  the  prior  pdf/a(o) .  This  result,  which 
follows  trivially  from  Eq.  (10),  was  first  given  by  Yao  [3], 

To  further  examine  the  consequences  of  the  estimator-correlator  structure,  we  now  assume  that 
T  >  0,  which  is  not  overly  restrictive  in  many  situations.  Immediately  we  may  observe  that  H0  is 
always  chosen  when  q0  <  q\.  Let  the  functions  g(q)  and  h(y)  be  defined  as 

s<?)  =  n^rW\  I  exp  H“}  Ua)da  ei) 


and 


h(y)  =  -  lny. 

With  these  definitions,  we  may  write 

F(q)  =  h(g(q))  +  C, 

where  F(q)  is  the  function  defined  in  Eq.  (14)  and  C  is  the  constant  h(g(0)} 
Straightforward  application  of  the  chain  rule  from  calculus  shows 


dq1 


dg_ 

dq 


d2g 

~8(q)-  8 


'd? 


g*(q) 


(22) 

(23) 


(24) 


Furthermore,  since  g(q )  is  positive,  we  may  show  that  this  second  derivative  is  <0  by  showing  that 
the  numerator  is  <0.  The  numerator  may  be  shown  to  be  <0  by  using  the  definition  of  the  function 
g(q)  given  by  Eq.  (21),  performing  the  indicated  differentiations,  and  applying  the  Schwartz  inequal¬ 
ity.  Thus  we  have  the  condition 


<  o.  (25) 

dq 
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At  this  point  we  assume  g0  <  q\  and  we  let /(g)  =  dF(q)/dq.  From  the  result  indicated  in  Eq.  (25),  it 
follows  that /(g)  is  a  monotonically  decreasing  function  of  q  and  that 

J<70 

M)dq-  (26) 

91 

Hence,  if  (go  -  gi)/(go)  >  T,  the  optimal  test  always  chooses  Hi  .  On  the  other  hand,  from  Eq.  (25)  it  also 
follows  that 

(*<70 

(go  -  gi)/(gi)  >  Kq)dq.  (27) 

Hence,  if  (g0  -  gi)/(gi)  <  T,  the  optimal  test  always  chooses  H0  .  From  these  considerations  we  are  led  to 
the  following  observation.  Whenever  the  detection  threshold  is  such  that  T  >  0,  then  the  optimal  detection 
procedure  may  be  implemented  in  a  sequential  manner  as  follows: 

1 .  If  go  <  gi ,  choose  Ho. 

2.  If  Eot**  l-^Kgo  —  gi)  T,  choose  H\. 

3.  If  E/fa  | x](g0  —  q\)  <  T,  choose  H0 . 

4.  Otherwise,  evaluate  the  log-likelihood  ratio  and  decide  accordingly. 

From  an  implementation  viewpoint,  these  considerations  may  reduct  the  computational  load  in  making 
a  decision,  since  the  full  log-likelihood  ratio  is  not  necessarily  evaluated  each  time  a  decision  is  reached. 

SUBOPTIMAL  DETECTION  STRUCTURES 

A  well-known  consequence  of  the  estimator-correlator  structure  is  its  implication  for  suboptimal 
detection.  If  the  optimal  estimator  E[a\z]  is  unavailable  or  difficult  to  work  with,  a  suboptimal  detector  can 
be  formulated  through  the  substitution  of  a  suboptimal  estimator  in  Eq.  (9).  For  example,  one  could  use  a 
parametric  representation  of  a  prior  density  with  a  sufficient  number  of  parameters  to  allow  a  large  number 
of  densities  to  be  fit  by  this  representation.  The  problem  is  then  reduced  to  a  problem  in  estimating  the 
parameters  of  this  prior  density  function,  a  task  that  could  be  implemented  adaptively. 

An  alternate  approach  to  suboptimal  detection  would  be  to  use  a  maximum  likelihood  estimator  (MLE) 
in  place  of  the  CME  in  the  estimator-correlator  structure.  From  the  standard  approach  to  finding  MLEs 
[37] ,  the  MLE  for  a  may  be  shown  to  be 


m 


“MLE - 

Z 


z 


=9i, 

2 


i  =  0.1. 


Use  of  this  estimator,  instead  of  the  CME,  in  Eq.  (9)  leads  to 


AmleW  = 


{x  -  5)'4>_1(JC  -  5) 


(28) 


(29) 
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If  we  expand  the  quadratic  form  in  the  denominator  and  divide  both  the  numerator  and  the  denominator  by 
the  numerator,  we  obtain 


AmleC*)  = 


2  Re(s'<f>  (x)  -  1s 


For  purposes  of  detection,  use  of  this  structure  is  equivalent  to  use  of  the  following  structure: 


(30) 


2  Re(sr$  -  sr< l>  '.v 

x‘<t>~'x 


(31) 


Except  that  our  result  is  for  complex  processes,  our  result  is  the  same  result  obtained  in  Ref.  7  by  an 
asymptotic  argument.  Therefore,  the  detector  that  is  shown  in  Ref.  7  to  be  asymptotically  optimal  is  shown 
here  to  be  given  by  a  suboptimal  implementation  of  the  estimator-correlator  with  an  MLE  substituted  for 
the  CME.  This  result  is  intuitively  satisfying  since  in  our  problem  we  expect  the  MLE  to  be  an  estimator  that 
is  asymptotically  equivalent  to  the  CME. 


A  problem  that  is  closely  related  to  ours  is  the  detection  of  signals  of  unknown  amplitude  and  phase  in 
Gaussian  noise  of  unknown  (but  nonrandom)  variance.  Korado  [14]  examined  this  problem  and  presented 
the  following  structure  as  optimal  for  this  problem: 


I  “I  I 

I  S  <P  AJC  | 

^xH^lx 


(32) 


If  we  take  Eq.  (29)  as  a  starting  point,  it  is  easy  to  show  that  the  MLE  of/3  (the  unknown  signal  amplitude  and 
phase)  is  given  by 


0MLE 


Substitution  of  this  result  into  Eq.  (29)  leads  to 


S  <P 


-1, 


5  $ 


-i: 


1- 


I  s  $ 


S 1  $  X '<f>  1 


For  purposes  of  detection,  this  structure  is  equivalent 


x  J 


(33) 


(34) 


Finally  a  rather  simplistic  estimator  is  to  estimate  a  by  a  constant  qq.  Use  of  this  estimator  in  the 
estimator-correlator  leads  to  the  following  test: 


A(.(a')  =  exp 


lx-(x—  s)'<I>  1  (jc  — 5-)] 


(35) 


which  easily  is  seen  to  lead  to  the  classical  matched  filter.  Clearly  if  the  pdf  of  a  is  <5(a  -  oq)  (i.e.,  the 
variance  is  a  known  constant),  then  the  estimator  is  correct  and  we  obtain  the  expected  result. 
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At  this  point  we  briefly  examine  the  structure  of  another  suboptimal  detector,  namely  the  locally 
optimal  detector  [15].  This  detector  is  given  in  general  by 


where  y  is  the  signal  amplitude.  By  straightforward  calculation,  one  may  show  that  for  our  problem 
this  structure  is  given  by 

A io(x)  =  £0[a|-*]  Re  (ej<t,s%~lx),  (37) 

i.e.,  the  locally  optimum  test  is  equal  to  the  normalized  matched  filter  multiplied  by  the  CME  of  a. 

Examination  of  the  suboptimal  tests  given  above  shows  them  all  to  be  essentially  matched  filters  that  are 
modified  in  some  way  to  account  for  the  unknown  or  statistical  behavior  of  a  .  This  observation  coupled  with 
our  result  that  the  optimal  detector  is  the  normalized  matched  filter  compared  to  an  adaptive  threshold  (i.e. , 
Eq.  (18))  leads  us  to  conjecture  that  good  suboptimal  tests  may  be  obtained  by  implementing  the  form 


Afo(*)  =  -T-  A(x) 

ay 


H i 

Refc'r1*)  >  fso(q0),  <38> 

Ho 

where/„,(i/o)  is  a  good  suboptimal  function  of  q  o-  However,  we  do  not  explore  this  issue  any  further  here. 

The  results  obtained  in  this  study  also  emphasize  a  point  about  the  detection  structures  for  these  types  of 
processes  (SIRPs).  Even  though  sample  vectors  from  these  processes  may  be  thought  of  as  samples  from  a 
Gaussian  process  with  an  unknown  and  varying  variance  [9],  the  detection  processing  should  involve 
estimation  of  the  quantity  a ,  not  of  the  quantity  r  (i .  e . ,  the  estimate  should  not  be  of  the  variance  directly) . 
Since  the  Gaussian  distribution  is  generally  parameterized  by  its  variance,  the  intuitive  notion  is  to  estimate 
the  u  lknown  variance;  the  results  here  show  that  this  intuitive  notion  does  not  lead  to  an  optimal  detection 
structure.  In  effect,  the  parameter  a  is  a  more  natural  parameter  than  the  variance  (a  conclusion  that  follows 
immediately  if  we  consider  the  Gaussian  distribution  in  the  context  of  the  exponential  family  with  the 
so-called  natural  parameterization). 

ALTERNATE  FORMULATION  OF  LIKELIHOOD  RATIO  AND  OPTIMAL  ESTIMATOR 

At  this  point  we  also  present  an  alternate  formulation  of  the  likelihood  ratio.  From  this  formulation,  we 
may  see  immediately  a  formulation  of  E[a  |  z]  that  does  not  require  explicit  knowledge  of  the  prior  pdf  fa(a) . 
Often,  a  prior  pdf,  being  a  pdf  of  a  parameter  rather  than  a  pdf  of  directly  observable  data,  is  difficult  to 
obtain.  Fortunately,  in  the  class  of  SIRPs,  an  alternate  pdf,  namely  the  marginal  amplitude  pdf  of  the 
complex  process,  may  be  used  to  solve  the  problem  without  knowledge  of  the  prior  pdf.  In  effect,  the 
marginal  amplitude  pdf  provides  information  that  is  equivalent  to  that  of  the  prior  pdf  for  this  class  of 
processes.  This  equivalence  is  fortunate,  since  in  practice  the  marginal  amplitude  pdf,  being  a  pdf  of  directly 
observable  data,  is  usually  easier  to  obtain  than  the  prior  pdf. 
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Since  we  are  examining  a  complex  SIRP,  we  have  also  under  hypothesis  H0  the  following  relationship: 


where 


*(M)  =  j0  oc\xj\  exp  fa(a)da,  j  =  1 


|  xj  |  is  the  amplitude  of  the  j th  component  of  x 
h(\xj  | )  =  pdf  of  \xj\,  j  =  \  ...  m. 


m 


(39) 


The  interpretation  of  Eq.  (39)  is  that  since  x  (under  Hq)  represents  a  vector  from  a  complex  zero- 
mean  Gaussian  random  process  with  random  variance,  the  marginal  pdf  of  the  amplitude  of  each  of 
the  components  of  this  vector  is  identical  and  Rayleigh  with  the  same  randomization  of  the  Rayleigh 
parameter  r  (=  1/a). 


From  Eq.  (39)  we  may  let  p  = 


and  rearrange  terms  to  obtain 


g(p )  =  =  jQ  exP  {~Pa}afo(pt)da. 


(40) 


Differentiating  with  respect  to  p ,  we  obtain  in  a  straightforward  fashion 


yfn- 1  p® 

C-1)"1'1  Up™71  =  Jo  exp  {~pa}da- 

Comparing  this  result  with  Eq.  (3),  we  obtain  finally 


(41) 


fi(x)  = 


(-D 


m-l 


dm 


i-l 


(27T)m|4>|  dp‘ 


m-l 


hd 2p) 

^2p 


q, 

P=Zi=- 


(42) 


Thus  we  see  that  knowledge  of  h(\xj\)  is  sufficient  to  represent  the  likelihood  ratio.  Often,  knowledge  of 
h(\xj\)  and  <f>  is  what  an  investigator  may  reasonably  assume;  for  instance  see  Ref.  9. 


To  show  how  this  result  relates  to  the  estimator-correlator  structure,  we  examine  the  CME  for  a: 

^  ”  (2xri*WM  I"  ^  exp  t  “ f} pWa-  (43) 

Note  that  Eq.  (43)  is  essentially  Eq.  3  in  which  we  have  substituted  m  +  1  for  m.  Therefore,  following  the 
same  derivation  that  leads  to  Eq.  (42)  yields 


pi  |  -j _  (~l)m  dm 

,Q!X  f(x)( :2*n$i  dPm 


KL 2p) 

'^2p 


p 


qr 

2 


(44) 
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Substitution  of  the  result  in  Eq.  (42)  into  Eq.  (44)  yields 


Ei\a  |x] 


f 


dm 

h(^2p) 

dpm 

^  J 

dm-\ 

h(^2p) 

dPm~X 

'G  J 

(45a) 


for  the  optimal  estimate  of  a  or 


E[a\p]  = 


dm 

h(^2p) 

dpm 

V2p 

dm~l 

"  h^2p) 

dpm~] 

^  J 

(45b) 


for  the  form  of  the  optimal  estimator  that  we  use  in  the  estimator-correlator.  Note  that  the  estimate  Ej[a  |jc]  is 
obtained  from  the  estimator  E[a\p]  by  evaluating  the  estimator  at  the  data  point  p  =  q-J 2. 

A  simple  check  shows  that  the  substitution  of  this  representation  of  E[a\p]  (i.e.,  £[a:|z])  into  Eq.  (9) 
leads  to  a  consistent  result.  As  we  stated  above,  this  result  is  important  because  in  practice  knowledge  of 
h(\xj\)  is  generally  available  whereas  knowledge  of fa(a)  is  generally  not  available.  However,  one  question 
that  remains  open  is  what  are  necessary  and  sufficient  conditions  on  /t(|jc,|)  such  that  it  has  the  representation 
given  in  Eq.  (39).  Such  knowledge  is  required  because  not  all  pdfs  h(\xj\)  have  such  a  representation.  Thus 
to  use  these  results  presently,  the  investigator  must  verify  in  some  way  that  the  h(\xj\)  of  interest  is 
compatible  with  such  a  representation. 

CONCLUSIONS 

In  this  report  we  reexamine  the  problem  of  signal  detection  in  the  class  of  spherically  invariant  random 
processes  and  show  that  this  problem  falls  within  the  framework  of  the  estimator-correlator  structure.  This 
structure  reveals  the  optimal  detector  to  be  a  function  of  the  optimal  estimator  of  a  random  quantity 
associated  with  the  spherically  invariant  random  process.  We  show  how  this  structure  leads  to  the 
interpretation  that  the  optimal  detector  is  an  adaptive  matched  filter.  We  then  demonstrate  how  some 
previously  obtained  results  are  easily  derived  as  suboptimal  implementations  of  the  estimator-correlator. 
Finally,  we  also  show  how  the  required  optimal  estimator  may  be  obtained  without  explicit  knowledge  of  a 
prior  pdf  provided  a  different  pdf,  namely  the  marginal  amplitude  pdf  of  the  complex  process,  which  is 
generally  easier  to  obtain  than  the  prior  pdf,  is  known. 

The  suboptimal  detectors  examined  here  may  be  classified  in  two  ways.  In  one  approach,  the  optimal 
detector  structure  is  retained  from  the  estimator-correlator,  but  a  suboptimal  estimator  is  used  in  place  of 
the  optimal  estimator.  In  the  other  approach,  the  optimal  estimator  is  retained  from  the  estimator- 
correlator,  but  the  function  of  this  estimator  that  forms  the  detector  is  different  from  that  of  the  optimal 
structure.  Interestingly,  each  of  these  suboptimal  approaches  results  in  a  detector  structure  that  is  in  essence 
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a  variation  on  the  classical  matched  filter.  Since  this  class  of  random  processes  is  a  generalization  of  the 
Gaussian  random  process,  this  result  is  intuitively  reasonable. 

Finally ,  since  this  class  of  processes  is  of  interest  as  a  model  for  some  types  of  practical  noise  processes, 
the  results  presented  here  should  lead  to  optimal  and  suboptimal  detection  schemes  for  practical  problems. 
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